A lattice Boltzmann study of reactive microflows 
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The role of geometrical micro-barriers on the conversion efficiency of reactive flows in narrow three- 
dimensional channels of millimetric size is investigated. Using a Lattice-Boltzmann-Lax-Wendroff 
code, we show that micro-barriers have an appreciable effect on the effective reaction efficiency of 
the device. If extrapolated to macroscopic scales, these effects can result in a sizeable increase of 
the overall reaction efficiency. 
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I. INTRODUCTION 



One of the outstanding frontiers of modern applied 
physics/mathematics consists in the formulation of mod- 
els and numerical tools for the description of complex 
phenomena involving multiple scales in space and time 
0. An important example of complex multiscale phe- 
nomena is the dynamics of reactive flows, a subject 
of wide interdisciplinary concern in theoretical and ap- 
plied science. The complexity of reactive flow dynamics 
is parametrized by three dimensionless quantities: the 
Reynolds number Re = UL/v, the Damkohler number 
Da — Th/Tc, and the Peclet number Pe = UH/D. Here 
U, L and H denote the macroscopic flow speed and lon- 
gitudinal/transversal lengths of the flow respectively, v is 
the fluid kinematic viscosity and D the pollutant molec- 
ular diffusivity. The quantities Tc and represent typi- 
cal timescales of chemical and hydrodynamic phenomena. 
High Reynolds numbers are associated with turbulence. 
High Damkohler numbers imply that chemistry is much 
faster than hydrodynamics, so that reactions are always 
in chemical equilibrium and take place in tiny regions. 
In the opposite regime the chemistry is slow and always 
takes place at local mechanical equilibrium. Finally, high 
Peclet numbers imply that the transported species stick 
tightly to the fluid carrier. Varying Re — Da — Pe and 
considering different device morphologies meets with an 
enormous variety of chemico-physical behaviours |Q. In 
this work we deal with low- Reynolds, fast-reacting flows 
with heterogeneus catalysis. In particular we wish to gain 
insights into the role of geometric micro-irregularities on 
the effective absorption rate of tracer species at catalytic 
boundaries. For a detailed study see also §. 



II. MATHEMATICAL MODEL OF REACTIVE 
MICROFLOW DYNAMICS 

We deal with a quasi-incompressible, isothermal flow 
with soluted species transported (advect and diffuse) by 



the flow and, upon reaching solid walls, they undergo 
catalytic chemical reactions. The basic equations of fluid 
motion are: 



(i) dtp + div pn — Q 
-VP+dw[2^(Vu4 



and (ii) dtpvL+div puu 
{Vuf)+\divu] 



where p is the flow density, u the flow speed, P — pT 
the fluid pressure, T the temperature and p,, A are the 
shear and bulk dynamic viscosities respectively (for the 
present case of quasi-incompressible flow with divM ~ 
the latter can safely be ignored). Finally, uu denotes the 
dyadic tensor UgU^,, a,b = x, y, z. 

Multispecies transport with chemical reactions is de- 
scribed by a generalized continuity-diffusion equation for 
each of s = 1, A^^ species: 



dtCs + divCgU ~ 

div[D,CTy{CJCT)] + (isSi^ - x^) , 



(1) 



where Cg denotes the mass density of the generic s-th 
species, Dg its mass diffusivity, Ct = J2s the total 
mass of transported species and fls is a chemical reaction 
term whose contribution is non-zero along the reactive 
surface described by the coordinate x^, {S{x) is the usual 
Dirac delta function). In the following the subscripts w 
and g mean "wall" (solid) and "gas" in a contact with 
the wall respectively. 

According to Fick's law, the outgoing (bulk-to-wall) 
diffusive mass flux (molecules per unit surface and time) 
is given by (hereafter species index s is omitted for sim- 
plicity): 



-Dd,C, 



±'~-'g\waU 



where d± denotes the normal-to-wall component of the 
gradient. Upon contact with solid walls, the transported 
species react according to the following empirical rate 
equation: 



dC^ 



Ja 



AS 

•av 



KrCu 



(2) 
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where AV^ is the volume element of the reactive wall and 
AS* is the surface element across which fluid-wall mass 
transfer takes place. In our case the ratio AV/AS is sim- 
ply the thickness of the reactive wall. Kc is the chem- 
ical reaction rate dictating species consumption once a 
molecule is absorbed by the wall. In the following we 
will use the common linear assumption 



AS* 



(3) 



where is the wall- fluid mass transfer rate. In prac- 
tice each boundary cell can be regarded as a microscopic 
chemical reactor sustained by the mass inflow from the 
fluid. Chemistry (Eqs. ||, ||) sets a time-scale for the 
steady-state mass exchange rate. At steady state we ob- 
tain: 



Hence 



A^ _ Cg 

•'^-^'"Ay-r^+r,' 

where = ^/K^, and Tc = 1/Kc- These expressions 
show that finite-rate chemistry [Kc > 0) ensures a non- 
zero steady wall outflux of pollutant. 




outlet 



inlet 



FIG. 1. Typical geometrical set up of the channel flow with 
a barrier on the bottom wall perpendicular to the flow of 
height h. 



As discussed in the introduction, they define the major 
dimensionless parameters 



Re = UH/v = ta/t^, Pe = UH/D = ta/td, 
Dac = Tc/TA, Da^^Tyj/TA- 



IV. CATALYTIC EFFICIENCY 



(4) 



III. THE COMPUTATIONAL METHOD 

The flow field is solved by a Lattice Boltzmann Equa- 
tion (LBE) method [§-0 while the multispecies transport 
and chemical reactions are handled with a variant of the 
Lax-Wendroff method (LW) ^ . The LW scheme repre- 
sents a numerically convenient choice recently developed 
to address multicomponent fluid transport (and reaction) 
within a LBE-like language. 

A. Multiscale considerations 

In this study, the unperturbed geometry of the cat- 
alytic device is a straight channel of size L lattice units 
along the flow direction (positive x direction) and H x H 
across it {y and z axes) . We add to this unperturbed ge- 
ometry a single protrusion (barrier) of unitary thickness 
at a fixed x — L/2 with height h in the z direction, and 
spanning the channel in the y direction (see Figj^). 

This problem involves at least four relevant time- 
scales. The relevant fluid scales are the advective and 
momentum-diffusive time: 

TA — L/U and t^, = H'^/v , 

The relevant time-scales for species dynamics are: 

TD = H^/D, T„ = K^^, and Tc = . 



The device efficiency is defined as amount of pollutant 
burned per unit mass injected: 

= — ' (5) 

where $(a;) = J[uC]{x, y, z)dydz is the longitudinal mass 
flow of the pollutant at section x and u is the x compo- 
nent of u [v and w will be the y and z components re- 
spectively). The in-out longitudinal flow deficit is equal 
to the amount T of pollutant absorbed at the catalytic 
wall per unit time. 

The goal of the optimization problem is to maximize 
r at a given <I>i„. This means maximizing complex 
configuration-dependent quantities, such as the wall dis- 
tribution of the pollutant and its normal-to-wall gradient. 
For future purposes, we find it convenient to recast the 
catalytic efficiency as 77 = 1 — T, where T is the channell 
transmittance T = ^out/^in- Roughly speaking, in the 
limit of fast-chemistry, this is controlled by the ratio of 
advection to diffusion timescales. It is intuitive that high 
efficiencies are associated with large values of the ratio 
ta/td, namely low-Peclet numbers. 

V. THE ROLE OF MICRO-IRREGULARITIES 

We now discuss the main qualitative effect of the 
micro-barrier from a microscopic point of view. 
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Firstly, it provides a potential enhancement of reac- 
tivity via the increase of the surface/volume ratio. How 
much of this potential is actually realized depends on the 
resulting flow configuration. 

Here, the fluid plays a two-faced role. First, geomet- 
rical restrictions lead to local fluid acceleration, hence 
less time for the pollutant to migrate from the bulk to 
the wall before being convected away by the mainstream 
flow. This effect may become appreciable on micro-scales 
for micro-flows with h/H ~ 0.1 (like in actual catalytic 
converters). Moreover, obstacles shield away part of the 
active surface (wake of the obstacle) where the fluid cir- 
culates at much reduced rates (stagnation) so that less 
pollutant is fed into the active surface. The size of the 
shielded region is proportional to the Reynolds number 
of the flow. On the other hand, if by some mechanism 
the flow proves capable of feeding the shielded region, 
then efficient absorption is restored simply because the 
pollutant is confined by recirculating patterns and has al- 
most infinite time to react without being convected away. 
This case is met mainly in presence of sufficiently ener- 
getic turbulent fluctuations at high values of the micro 
harrier-F eclei number Pen — ^jj^ 2> 1 where w' is the 
z-component of the velocity field at the barrier tip. 

With some appropriate approximations one can 
show that the efhcency is: 

r^o^l- e-^/' , (6) 

where I = llU/D, U = T^y^.^C / Y.y,.C , ll - 

CTH^/{2Cg) and r ~ {1/td + 1/(t, + T^)y\ 

Note that in the low absorbtion limit L <^ I, the above 
relation reduces to 770 — L/l, meaning that halving, say, 
the absorption length implies same efhciency with a twice 
shorter catalyzer. In the opposite high-absorption limit, 
L ^ I, the relative pay-off becomes increasingly less sig- 
nificant. 

We now turn to the case of a "perturbed" geometry. 
Let us begin by considering a single barrier of height h 
(Fig. [1]). The reference situation is a smooth channel at 
high Damkohler (Eq. ||). From we find an estimate 
of perturbative corrections in the smallness parameter 
g = h/H: 

^-Il^c^^^Re^,[Sc + K{a-l)] (7) 
770 2 H 

where A — H/L is the aspect ratio of the channel, Sc = 
v/D is the Schmidt number, and a is a regime dependent 
parameter. The wake length W can be estimated by 
W/h = KReu with K ~ 0.1. Three distinctive cases 
can be identified: (i) a = 0: the wake region is totally 
deactivated, absorption zero; (ii) a = 1: absorption in 
the wake region is exactly the same as for unperturbed 
flow; (iii) a > 1: the wake absorption is higher than 
with unperturbed ffow (back-flowing micro-vortices can 
hit the rear side of the barrier). 



VI. APPLICATION: REACTIVE FLOW OVER A 
MICROBARRIER 

The computational scheme has been applied to a fluid 
flowing in a millimeter-sized box of size 2x1x1 mil- 
limeters along the x, y, z directions with a perpendicular 
barrier of height h (see Fig. |l|). Upon using a 80 x 40 x 40 
grid, we obtain a lattice with dx = dy = dz = 0.0025 
{25 iim). We assume a real sound speed Vg = 300 m/s 
which becomes Cs = l/\/3 in lattice units. Therefore a 
time-step is equivalent to dt — Cg dx/Vg — 50ns. 

The flow is forced with a constant volumetric force 
which mimics the effects of a pressure gradient. The 
fluid flow carries a passive pollutant which is continu- 
ously injected at the inlet with a flat proflle across the 
channel. Diffusing across the flow, it reaches solid walls 
where it reacts according to a flrst order catalytic reac- 
tion: C + A P, where A denote an active catalyzer 
and P the reaction products. The initial conditions are: 

C{x, y, z) = 1 at the inlet, and = elsewhere 
p{x,y,z) ^ 1 

u{x, y, z) = Uo, v{x, y, z) = w{x, y,z)^0. 

The pollutant is then released at the open outlet, while 
flow periodicity is imposed at the inlet/outlet bound- 
aries. On the upper and lower walls, the flow speed is 
forced to vanish, whereas the fluid-wall mass exchange is 
modelled via a mass transfer rate equation of the form 
previously discussed. Our simulations refer to the fol- 
lowing values (in lattice units): Uo ~ 0.1 — 0.2, D — 0.1, 
v = 0.01, Kc = Kyj = 0.1. This implies Pe ~ 40, Re ~ 
400, Da > 80 (see also Eq. |). In order to study the 
effects of the barrier height h, we consider the follow- 
ing values: h = 0,2, 4, 8. The typical simulation time is 
t = 32000 time-steps (about 1.6 milliseconds in physical 
time) corresponding to two mass diffusion times across 
the channel. We may estimate the reference efficiency for 
the case of smooth channel: with U ~ 0.1, and t = 20, 
we obtain I ~ 200, hence 770 — 0.33. 

A typical two-dimensional cut of the flow pattern and 
pollutant spatial distribution in the section y — H/2 is 
shown in Figs. || and ||, which refer to the case h = 8. 

An extended (if feeble) recirculation pattern is well 
visible past the barrier. Also, enhanced concentration 
gradients on the tip of the barrier is easily recognized 
from Fig. ^. The integrated concentration of the pollu- 
tant C{x) = J2y z ^(^' ^) presented in Fig. ^ for the 
cases h = 0,2, 4, 8. The main highlight is a substantial 
reduction of the pollutant concentration with increasing 
barrier height. We measure also the the pollutant lon- 
gitudinal mass flow <I>(a;). The efficiency 77 is defined by 
Eq. |. The results are shown in Table ||, where subscript 
A refers to Eq. |^ with a = 1. These results are in a 
reasonable agreement with the analytical estimate apart 
deviations h = 8 for which the overall efficency is overes- 
timated. Leaving aside the initial portion of the channel. 
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u(x,z) at y=H/2:t=32000 
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FIG. 2. Typical two-dimensjpnal cut of the flow pattern 
with a single barrier of heigth h = S. Streamwise flow speed 
in the plane y = H/2. 
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0.295 


0.00 




R02 


1/20 


0.301 


0.02,0.025 




R04 


1/10 


0.312 


0.06,0.10 




R08 


2/10 


0.360 


0.22,0.40 
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TABLE I. Single barrier at x = 40: the effect of barrier 
20 Z height. 



C(x,z) at y=H/2: t=32000 
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FIG. 3. Concentration isocontjjurs with a single barrier of 
heigth /i = 8 on the plane y = _H'/2. 



Longitudinal pollutant concentration: single barrier 
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FIG. 4. Integrated concentration C(x) of the pollutant 
with a single barrier of height ft = 0, 2, 4, 8 after 32000 steps. 
The dashed line represent a theoretical evaluation with no 
barrier (ft = 0) and « ~ 200. 



streamwise speed U[Z) at x=3L/4,y=L/2 




U(Z) 

FIG. 5. Time evolution of the transversal streamwise speed 
u{z) at X = 3I//4 and y = L/2. Single barrier of varying 
zheight ft = 0,4,8 at f = 3200 and t = 32000. Note the 
backflow for ft = 8 at small z. 



our numerical data are pretty well fitted by an exponen- 
tial with absorption length I = 200, in a good agree- 
ment with the theoretical estimate I ~ 200. The barrier 
also promotes a potentially beneficial flow recirculation, 
which is well visible in Fig s. I and |. They clearly reveal 
a recirculating backflow for h = 8. For applications to 
many barriers see [||. 



VII. UPSCALING TO MACROSCOPIC DEVICES 

It is important to realize that even tiny improvements 
on the microscopic scale can result in pretty sizeable cu- 
mulative efi^ects on the macroscopic scale of the real de- 
vices, say 10 centimeters. The efficiency of an array of N 
serial micro-channels can be estimated simply as 



Vn 



N 



(8) 



It is readily recognized that even low single-channel ef- 
ficiencies can result in significant efficiencies of macro- 
scopic devices with = 10 — 100 (see Fig. 0). Equation 
^ with numerical data from present simulations provide 
satisfactory agreement with experimental data [p,p^. 

Nonetheless, extrapolations based on Eq. || must be 
taken very cautiously in the case of rugh or fractal walls 
[pil or of fully developed turbulence. 
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VIII. CONCLUSIONS 
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Although these simulations generally confirm qualita- 
tive expectations on the overall dependence on the ma- 
jor physical parameters, they also highlight the existence 
of non-pcrturbative effects, such as the onset of micro- 
vorticity in the wake of geometrical obstrusions, which 
are hardly amenable to analytical treatment. 

Work performed under NATO Grant PST.CLG. 
976357. 
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FIG. 6. Blow-up of the stres^pilines of the flow field past 
a barrier of height h = 8 located at a; = 40. The velocity 
direction in the closed streamlines of the vortex is clockwise. 
The recirculation effects are feeble and depletion is dominant. 
In fact for /i = 8 the local Peclet number is ~ 0.01-8/0.1 = 0.8, 
seemingly too small to provide micro-turbulent effects. 
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FIG. 7. Efficiency of a series of A'^ micro-channels as a func- 
tion of the single-channel transmittance. 
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